CERN-TH/2003-62 
DESY 03-029 



Reducing Residual-Mass Effects for Domain- Wall Fermions 

Karl Jansen a and Kei-ichi Nagai b 

& NIC/DESY Zeuthen, Platanenallee 6, D-15738 Zeuthen, Germany 
h CERN, Theory Division, CH-1211 Geneva 23, Switzerland 

Abstract 

It has been suggested to project out a number of low-lying eigenvalues of the four- 
dimensional Wilson-Dirac operator that generates the transfer matrix of domain-wall 
fermions in order to improve simulations with domain-wall fermions. We investigate how 
this projection method reduces the residual chiral symmetry-breaking effects for a Unite 
extent of the extra dimension. We use the standard Wilson as well as the renormalization- 
group-improved gauge action. In both cases we find a substantially reduced residual 
mass when the projection method is employed. In addition, the large fluctuations in this 
quantity disappear. 



1. Introduction 



Domain-wall fermions (DWF) preserve chiral symmetry ^ [21 13] when the lattice size in the 
5th direction, N s , is taken to infinity. The approach to the chiral limit is exponential in N s , 
with a rate given by the eigenvalues of the transfer matrix along the 5th direction, which is 
a local operator in 4 dimensions [3J 03 El EJ A measure of chiral symmetry breaking, taking 
place for finite N s , is the residual mass, m TCS , derived from the axial Ward-Takahashi identity. 

Even if the restoration of chiral symmetry is expected to be exponentially fast in N s , in 
practice m res can decrease very slowly as first shown by the CP-PACS collaboration [HUH]- 
The slow convergence of the residual mass is due to the existence of very small eigenvalues 
of the four-dimensional operator defining the transfer matrix along the 5th direction. In 
particular, at large N s these low-lying modes dominate the convergence rate [5] and render 
the recovery of chiral symmetry difficult. Even if the residual mass is very small, it is then not 
clear whether and what distortions of chiral symmetry are still present. Since large numerical 
simulations with DWF are being performed (see e.g. refs. [HIE! and the reviews I12j) it 
becomes important to find ways around this obstacle. Such solutions for improving the chiral 
properties of DWF then have to come from eliminating these low-lying modes. 

One idea to reduce these small eigenvalues is the improvement of the gauge actions [HI El 
I1U1 113j such as Iwasaki [Tl| or DBW2 ^21 ■ However, besides the potential difficulties with 
unitarity violations JB] an d the sampling problems of topological charge sectors J7j, this 
method does not solve the problem completely. For example, with the Iwasaki gauge action, 
the convergence rate also becomes slow at large N s The reason is that again very small 
eigenvalues of the transfer matrix appear in this case, though less frequently than for the 
Wilson gauge action. Using the DBW2 gauge action seems to be much better in this respect 
[T0| ITTj. but it is unclear whether these small eigenvalues could eventually appear there, too, 
leading to similar problems. A perturbative analysis JH] suggests a modification of the four- 
dimensional component of the domain wall operator to tackle the problem. This is, however, 
not yet tested in simulations. 

Another method to eliminate the disturbing effect of the small eigenvalues and the corre- 
sponding set of eigenstates of the transfer matrix is to project them out and lift them in a 
way that does not change the N s — > oo limit of the DWF operators |19l I2U| . In this paper, 
we investigate the projection method based on ref. JH1; where the projection is performed 
in the transfer matrix itself. In ref. |2Uj . an alternative projection is implemented through a 
modification of the boundary terms. The philosophy of both approaches is the same as the 
one using the transfer matrix. The aim of this article is to investigate the effects of the pro- 
jection method on the residual mass in quenched simulations. As we will see, the projection 
method works very well, leading to a substantial improvement in the residual mass. 

Let us emphasize that simulations with DWF can be considered under two aspects. The 
"purist's" approach demands exact chiral symmetry at non-zero lattice spacing. Here any 
violation of chiral symmetry (in practice up to machine precision) is not tolerable. Hence 
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the value of N s is to be taken as large as possible and the additive mass renormalization 
m rcs 7^ has to be eliminated. Thus the projection method discussed here, or any method 
leading to the same improvement, becomes an unavoidable necessity in this case. A different, 
more practical point of view is to consider DWF at finite, and even small, values of N s as a 
highly improved Wilson fermion. Also in this case, the projection method will accelerate the 
numerical simulation considerably and should therefore be employed. 

2. Domain- wall fermions and Ward-Takahashi identity 

In this section, we establish our notation and give the Ward-Takahashi identity in order to 
define the residual mass. For completeness, we give here the definition of the domain-wall 
operator and its relation to the 4D operator satisfying the Ginsparg- Wilson equation |21j . 
We follow the presentation of ^Hj . Derivations of this formulae can be found in [21 01 HI 13 Hj • 
The 5D domain-wall operator is defined as 

T> = ^H(d* s +d s )-a s d* s d s } + M , (2.1) 

where s denotes a lattice site in the 5th direction (1 < s < N s ), a s is the corresponding lattice 
spacing, and d* and d s are the free forward and backward derivatives. 

The operator M is obtained from the standard 4D Wilson-Dirac operator by 

M = D w - m (2.2) 

with 

Ay = ^{>( V A* + V m) ~ a K V »} . (2.3) 

Here and are the gauge covariant forward and backward derivatives and a is the lattice 
spacing in the four physical dimensions /i = 1, . . . , 4. The domain- wall parameter mo obeys 

< a s m < 2 , < am < 2 . (2.4) 

Note that the lattice spacings a s and a can be different in general. The boundary conditions 
in the DWF formulation in the 5th direction is 

P+V>(0, x) = P-^{N S + 1, x) = , (2.5) 

where P± = |(1 ± 75). In these settings, the chiral modes with opposite chiralities are 
localized on 4D boundary planes at s = 1 and s = N s . 

The 4D quark fields are constructed from the left and right boundary (chiral) modes, as 
follows: 

q(x) = P^(l,x)+P + ^(N s ,x) , q(x) = ij(l,x)P + +ij(N s ,x)P_ . (2.6) 
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A bare quark mass term is introduced by adding to eq. (|2.1|) the term 

m f {4>(l,x)P + Tp(N s ,x) +$(N s ,x)P-'i(j(1,x)} = m f q(x)q(x) . (2.7) 
The propagator of the quark fields is related to an effective 4D operator D]v g [3JIH1II1 

(q(x)q(x)) = 2 -^^ , (2.8) 

with 

DN s ,m f = (1 - am/) D Ns + 2m/ . (2.9) 
In terms of the operators K± , 

,,11 a s M 

K ± = -±-ys — 1 — r-7 , 2.10 
± 2 2 lb 2 + a s .M ' v ; 

Z?at s is given by 

K N S _ K N S 

aD ^= l + ^ K + N s+K N B ■ (2-11) 

From this equation, it is easy to show that 

aD = lim aD Ns = 1+75 sign (K+ - KJ) , (2.12) 

which is written as 

00 = ' (2 - 13) 

A = - aM KA . (2.14) 
2 + a s ,M v ' 



The operator D in eq. (|2.13|) satisfies the Ginsparg-Wilson relation. The only difference to 
Neuberger's operator |22j is the definition of A. Neuberger's operator is obtained from eqs. 



(|2.13j) and 1)2. 14|) by taking the limit a s — > 0. 

In the limit N s — > 00, the 5D formulation of DWF is completely equivalent to a 4D lattice 
formulation of Ginsparg-Wilson fermions satisfying an exact chiral symmetry. However, in 
a realistic simulation N s is kept finite, of course. In this situation, the chiral symmetry is 
explicitly broken by the residual terms 5D = Dn s — D. A measurement of the effects of 
this chiral symmetry breaking is the so-called residual mass m TCS , derived from the axial 
Ward-Takahashi identity. The chiral transformation of DWF is defined as 

Sip(s,x) = iQ(s)eip(s,x) , 5ijj(s,x) = —i'iJj(s,x)Q(s)e, (2.15) 

where Q(s) = sign(A^ s — 2s + 1) and e is an infinitesimal transformation parameter. Under 
this transformation, the quark fields are transformed as the usual chiral transformation: 
5q(x) = i"f^eq{x) , 5q(x) = iq{x)^§e . Therefore the axial Ward-Takahashi identity is 

X;(V M A M (x)P(0)) = 2m f (P(x)P(0)) + 2(J 5g (x)P(0)) , (2.16) 
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where A fM (x) is the axial- vector current and P(x) is the pseudo-scalar density; A )M (x) and 
P(x) are given as 

A IM (x)=J2Q(s)Ms,x) , P(x) = q{x) lb q{x) , (2.17) 

s 

where 

M x ) = 2 [^ s ' " T^WC*. x + »)~ fa, x + M )(l + 7 M )£/t(x)^(s, x)] . (2.18) 
The additional term Js g represents the explicit breaking of chiral symmetry, 

J 59 = -^(y> x ) p -^(^ + 1 > x ) +$(y + 1 ' x ) p+ ^(t> x ) ■ (2,19) 

For a smooth gauge field background, the term (J^qP) vanishes exponentially fast [3J HE HS1 
as N s is increased. In realistic simulations, however, the gauge fields can be rough and it may 
happen that the rate of convergence in N s is rather poor. The breaking of chiral symmetry 
can be quantified by the values of m res - Let us define the ratio R(t): 

_ Ex(^(x,t)P(o,Q)) f2 . 

m ~ L(P(x,i)P(o,o)) • (2 ' 20) 

The usual definition of m res is the average of the quantity R(t) at large time separations. A 
necessary but, maybe, not sufficient condition to recover fully on-shell chiral symmetry at a 
non-vanishing value of the lattice spacing is that this quantity be negligible. 



3. Eigenvalues of A^A 

For gauge configurations with a restricted value of the plaquette (so-called admissible config- 
urations) [SJ[23], the operator A^A has been shown to have a spectral gap, < u < A^A < v, 
ensuring the exponential suppression of the residual mass in N s . However, in realistic simu- 
lations as performed today, the plaquette bound is not satisfied and it is important to study 
the distribution of the eigenvalues of A^A in numerical simulations. 

The eigenvalues of A^A can be obtained through the generalized 4D eigenvalue equation 

m 

agM^Mip = A (2 + a s M)^(2 + a s M)tp . (3.1) 

The low-lying (maximal) eigenvalues can be computed by minimizing (maximizing) the gen- 
eralized Ritz functional 1 

(^\(2 + a s My(2 + a s MM 1 ' ' 

The interested reader may obtain more details on request. 
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using a straightforward generalization of the algorithm described in ref. [22]. Notice that in 
this method no inversion of the matrix (2 + a s A4)^(2 + a s A4) is needed. Higher eigenvalues 
can be calculated by modifying the operator A4'A4 in the numerator in eq. (j3.2j) . so that 
the already computed eigenvalues are shifted to larger values. This can be achieved |19j by 
substituting 

M'M M^M + M^M\A)(A\(2 + a s M)\2 + a s M) , (3.3) 

Xi, ipi being the already computed (lower) eigenvalues and eigenvectors. 

Figure^shows the eleven lowest eigenvalues of the operator A^A as a function of the Monte 
Carlo time (*mc)- Here and throughout the paper we use the quenched approximation and 
set a s = 1. The data in Fig. ^are obtained with the Wilson gauge action at (3 = 6.0 on a 
12 3 x 24 lattice, setting tuq = 1.8. As expected, very small eigenvalues appear frequently. 
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Figure 1: The 11 lowest eigenvalues of the operator A^A as a function of Monte Carlo time 
iMC at /? = 6.0 and mo = 1.8 on a 12 3 x 24 lattice. The open diamonds denote the lowest 
eigenvalue. 



The minimum rate of convergence in N s of the operator D^ 3 is given by 



ijj = mm \Uj 



Ui = In ■ 



(3.4) 



where \ are the eigenvalues of A^A [Ej. Figure El shows the inverse convergence rate com- 
puted from the eigenvalues in Fig. ^ Clearly, the low-lying eigenvalues of A^A lead to a slow 
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convergence, causing the simulation to become expensive. 

10 3 , 




Figure 2: The inverse convergence rate u i 1 computed from the eigenvalues plotted in Fig^ 

We also explored the eigenvalues for other gauge actions such as the Iwasaki |14j and 
DBW2 |15j ones. An example for these eigenvalues is plotted in Fig. |31 In that figure we 
average over 20 gauge configurations. The parameters of the gauge actions were chosen such 
that in each case the value of the lattice spacing is a = 0.093 fm, leading to setting f3 = 6.0 
for the Wilson action, (3 = 2.6 for the Iwasaki one and (5 = 1.04 for the DBW2 one. Since 
also the lattice size was fixed to be 12 3 x 24 we have for the different gauge actions the same 
physical situation. For the Wilson action we observe small values for the lowest-lying modes. 
This is improved substantially by employing the Iwasaki action and even more when using 
the DBW2 action. Note that the 11th low- lying eigenvalue of the Wilson action corresponds 
to the lowest eigenvalue of the Iwasaki action. We checked for the Wilson and the Iwasaki 
action that this picture does not change when we decrease the value of the lattice spacing 
down to a = 0.05 fm. This confirms that the convergence in N s is faster when the gauge 
action is improved [HI 1101 HHj . As we will see below a conclusion that improved gauge actions 
by themselves would completely cure the problem of a slow convergence rate is premature, 
however. 
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Figure 3: Averaged eigenvalues for the Wilson, Iwasaki and DBW2 gauge actions as a function 
of the eigenvalue number. The lattice spacing a = 0.093 fm used is the same for all gauge 
actions. 

4. Improvement of domain-wall fermion 

The decay rate of the residual mass in N s is controlled by the small eigenvalues of A* A. 
For the Wilson gauge action very small eigenvalues occur, leading to a slow convergence. 
Although the situation is improved for the Iwasaki gauge action, as we saw above, it was 
observed that even in this case for large values of N s the convergence turned to become very 
slow [B1IH1- It thus seems to be necessary to test methods as proposed in JH] that modify 
the fermionic part of the DWF action by projecting out the small eigenvalues of A' A. These 
methods can be used alternatively -or even in addition- to employing improved gauge actions. 
The key observation in JS] is that the relations in eqs. (|2,13|) and (|2.14|l hold true for any 
choice of M as long as 

M ] = 75M75 , det(2 + a s M) / . (4.1) 

This fact may be used to construct an improved M for which the very low eigenvalues of 
A' A disappear. 
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Let us, for completeness, repeat the construction of the improved operator here again 
following ^5] . The basic idea is to find the new operator Ai satisfying the following relation; 

~ = *sM = A + V (Sfc - a k )^v k ® v{ , (4.2) 

2 + a s M ^ 

where v k is the eigenvector of the following equation 

-f 5 Av k = a k v k , k = 1, . . . , r , (v k ,vi) = 5 k[ . (4.3) 

Therefore an improved DWF operator, D^%, can be obtained from eq. (|2.1j) after substituting 
M with M defined as 

r 

a s M = a s M - ^ X kl w k <g> w}j 5 , (4.4) 

k,l=l 

where 

w k = {2 + a s M)j 5 v k (4.5) 

and 

{X~ l ) kl = 25 kl {a k -a k )- 1 + (v k , Wl ) . (4.6) 

It is easy to see that 75^4 has the same eigenvectors as 75 A; however, all eigenvalues a k , 
k = 1, . . . , r, are replaced by a k . The limit N s — > 00 of Dn s is of course unchanged by this 
modification, provided 

sign(2 fc ) = sign(a fc ) . (4.7) 
The choice of \a k \ is not unique. We will choose here 

a k = 2 sign(a fc )|a ; | , 1 < k < r = /c max , (4.8) 

where A; max is the number of eigenvalues projected out and I can be chosen freely. A natural 
choice is / = k mSuX such that all low-lying eigenvalues are moved to be twice higher than the 
largest eigenvalue projected out. We also tried, however, different values of I and found that 
the improvement is not very sensitive to the precise choice of a k , provided it is larger than 

OLh 

<Tnax 

Our statistics is typically 60 configurations for the Wilson gauge action and 20 config- 
urations for the Iwasaki action. We did not explore the DBW2 action extensively. The 
parameters of the gauge actions were chosen as before such that o" 1 = 2 GeV, which means 
a choice of (3 = 6.0 for the Wilson gauge action and (3 = 2.6 for the Iwasaki one. The lattice 
sizes were 12 3 x 24 x N s and 16 3 x 24 x N s for the two actions, respectively. The domain-wall 
mass was tuq = 1.8 and we worked at a quark mass of mt = 0.02. 

We have measured the residual mass from R{t) in eq. (|2.20|) as the average of R{t) for t 
typically in the interval 4 < t < 20 for a time extent of the lattice of T = 24. R(t) is shown 
in Fig. for the case when no projection is performed. For each value of N s we have the 
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Figure 4: The ratio R(t) = (ppj as a function of Euclidean time. The Wilson gauge action 
is chosen with simulation parameters as given in the text. No projection of eigenvalues is 
performed. 



same statistics. Although, with increasing N s , the residual mass m res decreases, it does so 
rather slowly; furthermore, as N s increases, large fluctuations in R(t) occur, rendering the 
determination of the residual mass difficult. These large fluctuations also suggest that the 
residual chirality-breaking effects in other quantities might be very hard to estimate, taking 
only m res as a measure of these effects. 

In Fig. we show R(t) when we project out a number of eigenvalues. As expected, the 
projection of the low eigenvalues decreases the residual mass significantly with respect to 
Fig.H 

The more eigenvalues are projected the smaller the residual mass is. Another important 
feature is that the fluctuations in R(t) become much smaller when a sufficiently large number 
of eigenvalues is projected out; in this case 10 seems to be a good choice. This is very clearly 
seen in Fig. El where we show the value of the ratio R(t), 

E x P(x,t)P(o,o) 

computed on single configurations at t = 12 as a function of Monte Carlo time. The spikes 
are substantially damped with the projection. Finally, when the projection is implemented, 
the decrease of the residual mass with N s is much faster. 
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Figure 5: Same as Fig. |I] but now with the projection of eigenvalues employed. 
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Figure 6: The quantity R(t) — (R(t)) at t = 12, see eq. (|4.9j) . as a function of Monte Carlo 
time for the Wilson gauge action, with and without projection. 
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In summary, it is clear that the projection method has a drastic effect on the value and dis- 
persion of the residual mass. However a sufficient (O(10) in our setup) number of eigenvalues 
have to be projected out. 



5. DWF with improved gauge actions 

In our simulations with the Iwasaki gauge action we found, even in our small sample of only 
20 configurations, very low-lying eigenvalues of A' A. In order to see the effect of these modes 
we plot, in Fig. [7J the ratio R(t) of eq. (|4.9j) . for two of these configurations (note that no 
averaging is involved here). The figures indicate that we will find, also for the Iwasaki gauge 
action, the same problem as for the Wilson gauge action. When no projection is performed, 
the correlation function shows a spiky behaviour, which may lead to large fluctuations in 
R(t) and hence to a very difficult determination of the residual mass. This is also confirmed 
in the ratio of averaged values, R(t), as shown in Fig. (HI The pattern resembles the case of 
the Wilson gauge action. For small values of N s the effect of the projection is not noticeable. 
For larger values of N s , we see that R(t) is lowered when the eigenvalues are projected out 
and that the fluctuations of this quantity are strongly damped. 
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Figure 7: R(t) as a function of time in the Iwasaki action. The circles show the results for 
the no-projection case and the diamonds for the case when 3 eigenvalues are projected out. 



We also made an attempt to see how the projection method affects R(t) for the DBW2 
action. For the simulations we chose j3 = 1.04, which corresponds again to a" 1 = 2 GeV. 
Thus we study the same physical situation with the Wilson and Iwasaki gauge actions. The 
lattice size was chosen to be 16 3 x 32 x N s and mo = 1.7. The fermion mass was taken to be 
m f = 0.02. 

We observe in Fig. |§]that the residual mass is not changed very much by the projection. 
We attribute this to the fact that in our small statistical sample no very low-lying eigenvalues 
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Figure 8: T&jff without and with projection for the Iwasaki gauge action. 



of A^A could be detected. We see, however, from the same figure that the statistical error is 
substantially reduced for certain values of t when the projection of eigenvalues is employed. 

The fact that R(t) shows large fluctuations, even though there are no very small low-lying 
eigenvalues, points toward the suspicion that also the eigenvectors may play an important role. 
In particular the localization properties of these eigenmodes may lead to large fluctuations 
as discussed in j^B]. Although this point deserves further investigation, we did not perform 
such a study here. To conclude, from a negligible average value of the residual mass, that 
chiral symmetry is restored is certainly questionable when the dispersion of the residual mass 
is large and not gaussian. A much safer situation would be to ensure that the residual mass 
is bounded from above for all configurations. The projection method ensures that this is the 
case. 

To summarize, in Fig. ^H we show the comparison of the behaviour of the residual mass as 
a function of N s for different gauge actions and for different numbers of projected eigenvalues. 
For a fixed gauge action, we find that at small N s there is almost no effect from the projection 
method. 

This can be explained by a simple qualitative argument with the formula suggested in 

OB 03 mi; 



-a k N s 



dap(a)e~ aNs , (5.1) 
where p(a) is the eigenvalue density in the continuum. This (qualitative) formula describes 
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Figure 9: ^ipp\ without and with projection at N s 



16 and 32, for the DBW2 gauge action. 



the behaviour of function of N s . The formula contains two factors, the eigenvalue 

density and the exponential supression factor e~ aNs . For small values of N s , not only do the 
low-lying modes contribute to the sum in eq. ()5.1j) . but also the bulk modes since they are 
not supressed sufficiently. When projecting out a few number of low-lying eigenmodes, the 
eigenvalue density and the exponential factor remain almost unchanged and hence also the 
residual mass is not affected very much for small values of N s . In such a case, it would be 
necessary to project out a large number of eigenmodes to make m res decrease. When N s is 
chosen to be large, on the other hand, the bulk mode contributions to the sum in eq. (|5.1|) 
will die out and only the small eigenvalue contributions will survive. As a consequence, the 
factor e~ aNs becomes much smaller after projecting out even only a few low-lying (isolated) 
eigenmodes. This should hence lead to a large improvement, i.e. a substantial decrease of 
the residual mass when the projection method is active. As Fig. EH clearly shows, this is 
indeed the case. For the Wilson gauge action at N s = 48, the value of the residual mass is 
decreased by several orders of magnitude when 10 eigenvalues are projected out. We made 
a rough check for the Iwasaki gauge action that also in this case the residual mass decreases 
substantially, choosing N s = 40. Thus the very slow decrease of the residual mass as a 
function of N s in the original DWF formulation with no projection is cured by projecting out 
a few O(10) eigenvalues. 
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Figure 10: A compilation of the residual mass as a function of N s for various gauge actions 
and various choices of projecting eigenvalues. The filled symbols correspond to our own 
results. The data of the DBW2 action are taken from ^U] and the ones for the Iwasaki action 
from [SI The lines are just to guide the eye. 



6. Conclusion 



We have studied the effect of modifying the fermion action of DWF by projecting out a few 
low- lying eigenvalues of the underlying transfer matrix ^5]. By measuring the correlation 
function leading to a determination of the residual mass and the residual mass itself as a 
function of N s , we find a significant improvement in the restoration of chiral symmetry for 
quenched DWF at large N s . 

The reason is that in the large- N s limit the low-lying eigenvalues of A^A are responsible for 
the exponential convergence rate of DWF in N s to its chiral invariant limit. These eigenvalues 
then dominate the behaviour of the residual mass and whenever very small low-lying modes 
appear they lead to a very slow decrease of the residual mass as N s is increased. Projecting 
out a small number of these modes can therefore help considerably to lower the values of the 
residual mass. We have confirmed this picture in practical simulations, using the Wilson and 
the Iwasaki gauge actions. We observe that when a sufficient number, i.e. O(10), eigenvalues 
are projected out, the residual mass vanishes rapidly with increasing N s . 

Let us end our discussion with three remarks. 
(i) Projecting out a number of low-lying eigenvalues shows a strong effect not only on the 
value but also on the fluctuations of the correlation function R(t) in eq. 1)2. 20j) and hence 
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of the residual mass. The damping of the fluctuations takes place even when no very small 
eigenvalues occur in the simulation, as in the case of the DBW2 action. It thus seems that 
also the eigenvectors and in particular their localization properties play an important role. It 
is unclear to us, and we did not investigate this here, how far also other correlation functions 
are affected by this phenomenon. One possible explanation [22] relies on the relation of the 
eigenvalues and eigenvectors of D\y and D^ ajrrlf . The study of this correspondence clearly 
deserves further efforts using non-perturbative methods. 

(ii) The method of projecting out eigenvalues as studied here can be used on top of other 
improvements such as using improved gauge actions or improved fermion actions. The pro- 
jection method is not very costly and produces only a small numerical overhead. Thus we 
advocate to employ the projection method in any simulation done with DWF. 
(Hi) We expect that the projection of the low-lying eigenvalues should play an even more 
important role in the case of dynamical simulations with DWF as the behaviour of the 5D 
fermionic kernel will be affected by the problems discussed in (i), too. We envisage that 
such a dynamical computation with the projection of low-lying eigenvalues can be performed 
along the lines of refs. |27| 1281 1291 I30| by estimating the full DWF operator stochastically. In 
this case the projection can be done easily. 
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